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1. Introduction and stability result 

1.1. A spatial ALOHA 

We consider a random spatial service system governed by an ALOHA-type algo- 
rithm. More precisely, time is slotted, during a time slot n — 1,2,.. a random number 

of users arrive in the system, and at each slot every user in the system requires 
service (transmission) with a certain (transmission) probability independently of all 
other users. The sequence is assumed to be i.i.d. 

The ALOHA multi-access algorithm was first proposed by Abramson lj. The 
slotted scheme was introduced by Roberts [15] ■ We consider the latter setting. In 
the conventional centralised slotted ALOHA model, there is a single server. At the 
beginning of a time slot n, the total number W n of users in the system is known, and 
each of them asks for service (transmission) with probability independently of all 
others. If two or more users attempt transmissions simultaneously, then transmissions 
collide, the users remain in the system and try to transmit later. All service times are 
equal to 1, and the server is free at the beginning of any time slot. It is easy to show 
that for this system the maximum throughput is equal to e . Further, the Markov 
chain is positive recurrent if E£i < e _1 and transient if E£i > e _1 . 

When the information on the numbers W n of users is unavailable, various decen- 
tralised adaptive algorithms have been introduced and studied. Algorithms of this type 
use the information on what has occurred in the previous time slot: either a conflict or 
a successful transmission or no attempted transmissions. More precisely, let B n be the 
number of users trying to transmit during a time slot n. In a decentralised algorithm, 
only the values of min{i?„, 2} are available at time n+1. Assuming that at any time slot 
n the transmission probability is the same for all users, Hajek jl lj proved that if {£ ra } 
are i.i.d. with a finite exponential moment, then E£ n < e _1 ~ 0.37 is still necessary 
and sufficient for the existence of a stable algorithm. Mikhailov [T3] generalised this 
result by weakening the exponential moment assumption to the requirement for only 
the second moment to exist, while Foss [5] generalised it further by dropping this as 
well as the independence assumption. We also refer to Ephremides and Hajek [6J for a 
survey which includes, in particular, results in this direction. 

The model described above ignores the network's spatial diversity. In this paper we 
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present a new model which allows the possibility for interaction between users to occur 
only when the distance between these users is small. A limited spatial interaction 
is a common feature of wireless networks and is thus very important for practical 
consideration. 

In this paper we consider only spatial centralised schemes (protocols) where the total 
number of users in a neighbourhood is known. The study of decentralised protocols is 
a subject of our future research. Our main result establishes the stability of the model 
under consideration. For that, we use the fluid approximation approach. We show that 
the state is repelled by the boundary, and in the interior, the Euclidean norm turns 
out to be a Lyapunov function for the fluid model. Other results of the paper include 
convergence rates to stationarity and local (in)stability. 

The remainder of the paper is organised as follows. The rest of this introduction is 
devoted to the description of the model and the statement of our main stability result. 
In Section [5] we prove that fluid limits of the workload in the system satisfy a certain 
differential equation. Section [3] is devoted to the analysis of the behaviour of fluid 
limits on the boundary of the positive orthant. In Section 0] we present a proof of our 
main stability result and formulate one of its possible generalisations. Section [5] deals 
with rates of convergence towards the steady state. Sections [6] and [7] contain some 
interesting results on the behaviour of solutions to the differential equation. Finally, in 
Section [8] we conclude with some extensions of our model, which arc in a certain sense 
more applicable to real systems. In particular, these extensions include systems where 
various changes in environment conditions may result in changes in the radius and/or 
direction of interference between the message transmissions). 

1.2. Model description 

Let Q = (V, £) be an undirected graph with a finite set of vertices, say V = {1, K}. 
We suppose that Q is connected. For the graph Q we use the standard notion of the 
graph distance. For i £ V, let V% — {i} U {j £ V : such that £ £} be the 

neighbourhood of a vertex i, i.e. the set of vertices at a maximum distance of 1 from 
i in the graph. In the particular case where all vertices have the same number of 
neighbours, we denote by V the cardinality of Vi, card(Vi) = V. 

We introduce now a service system with spatial (neighbourhood) interactions as- 
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sociated with the graph Q. We assume that time is slotted, i.e., arrivals and services 
may occur only at times n = 1,2,.... Suppose that there is a service station at 
each vertex of Q. The arrival process is denoted by A = (A(n)) ne ^, where A(n) = 
(Ai(n), Ajf(n)) e N K , and Ai(n) is the number of users arriving at time n at a 



number of users arriving between time instances s and t— . We suppose that (A(n)) is 
an i.i.d. sequence (while the coordinates of the vectors may be dependent). We also 
suppose that EAj(n) = Xi > for i = 1, . . . , K. If values of \ do not depend on i, 
then we write Xi = A, i = 1, . . . , K. 

Let W{n) — (Wi(n), Wk(ti)) E R+ , where Wi(n) is the number of users at vertex 
i at time n. At time n, a user at vertex i attempts transmission independently of the 
others with probability 1 / ^ e v< ( n ) • This user receives service if he is the only user 
attempting transmission in Vi at time n. We suppose that all service times are equal to 
1 and that any user leaves the system immediately upon service completion. Let N t (n) 
be the number of users attempting transmission at time n at vertex i. Then iV,(n) is 

a binomial random variable with parameters [ WAri), — — ) and {NAn), 1 < 

V 22jev t w i( n ) J 

i < K} are independent variables conditioned on W(n). The sequence {VF(ri)} forms 

a time-homogeneous irreducible Markov Chain, which satisfies the recursion 



W i {n) = W i (n-l) + A i {n)-l(N i (n-l) = l) ] | l{N j {n - 1) = 0). (1) 



To explicitly show the dependence of W{n) on the initial condition W(0) — x, we 
may sometimes write W x (n). 

If Xi > 0, then the i-th component of the drift vector is given by the expression 



vertex i. For s < t, denote by A(s, t) = J^, 



■ rtl-i 

n= \s~\ 



A{n) the vector representing the 



E [Wi(l)-Wi(0)\W(0) = x] = Xi 




and if x t = 0, then E [ W t {\) - Wi(0)\W{0) = x ]= \. 



We rewrite the expression for the drift vector in the following way: 



E [W(l) - W(0)|W(0) = x]= A- G(x). 
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Here A is a X-dimcnsional vector with its z-th component being equal to A, and G(x) = 
(Gi(x), . . . , Gk{x)) is a function from R K to H K defined by 

qj x \ = { ** ( X ~ E*^**) IW«\{i} ( x - t^l,**) ' if Xi > °' 
[o, if Xj=0. 

For a; G E 1 ^, we denote 0i(a;) = ^ Xi — . Let = (0i(x), ^k(x)). Note that G,: 
is bounded from above by 1 and if X^fcev ^fc > then 

lim Gi{tx) = Gi(x) = <f>i(x)e~ 

t — ' + 00 

In particular, Gi is a homogeneous function of order 0, i.e. Gi{cx) = Gi(x) for any 
c> 0. 

We now comment on the model. In this paper, we mostly consider the symmetric 
case, where A^ = A for all i = 1,...,K and the graph Q is (V — l)-regular: the 
cardinality of Vi is equal to V for all i. Notice that, even in this case, the graph Q is 
not necessarily completely symmetric. 

Note also that the system is not monotone. Indeed, x < y (component-wise) does 
not imply that W v (l) stochastically dominates W x (l). Neither is the system monotone 
with respect to the graph structure: if G\ is embedded into G*2, this does not imply that 
the workload process built on graph G\ is stochastically dominated by the workload 
built on graph G^. 

We also present a number of generalisations to a non-symmetric case. In partic- 
ular, using methods suggested recently in [18] , we formulate in Remark |3] sufficient 
conditions for the stability of the system with space-inhomogeneous input. Some other 
generalisations of the model are proposed in Section [51 

1.3. Stability result 

We first explain the intuition hidden behind the result. Consider the symmetric 
case, Aj = A, card{Vi) = V, i = 1, . . . , K . 

The access protocol favours an equilibrium of the workload in the network: assume 
that the workload at node i is much larger than the workload in its neighbouring 
nodes, Vi \ {i}. Then <f>i(x) is close to 1, whereas (f>j(x) is close to for all the other 
nodes j in Vi. Thus the workload at node j in Vi will tend to get closer to the 
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workload at node i. From this balance mechanism one can guess that the diagonal 
A = {.i e M. K : x\ — X2--- — xk} is an attractive set. 

If the workload vector belongs to the diagonal: W(0) = cl where c € N, then we 
obtain: 



Hence, as c tends to infinity, the drift vector converges to (A — e 1 /V)l. 

So finally, we end up with the conjecture that if A < e _1 /V, then the Markov chain 
W is ergodic. This conjecture is clearly true for the fully isolated graph and for the 
complete graph. 

The reasons that led to this conjecture appear to be wrong (this will follow from 
the results of Sections [5] and which show that in general the main diagonal is not 
attractive). However, the conjecture itself is true and we can formulate our main 
stability result that will be proved in Section |U 

Theorem 1. Assume the model to be symmetric. 

(i) If A < e^ 1 /V, then the Markov chain W is positive recurrent. 

(ii) Assume further that P(^i(l) = A 2 (l) = ... = A K (l) = 0) > 0. Then W is 
ergodic, i. e. there exists a unique stationary distribution and the distribution of W n 
converges to the stationary one in the total variation norm. 

Our proof of (i) is based on the fluid approximation approach. We will show that 
all fluid limits satisfy a certain differential equation and then proceed with the analysis 
of that equation. The proof of (ii) is standard: the state = (0, 0, 0) is achievable 
from any other state and from itself. So, the state is positive recurrent (due to (i)) 
and the Markov chain W is aperiodic. The ergodicity follows. 

Our heuristics suggest also that if A > e~ 1 /V then W is transient. Corollary |4] in 
Section [6] is a partial result which corroborates this intuition. 




2. Fluid approximation method 



This section deals with a general (not necessarily symmetric) graph. 
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2.1. General properties 

In what follows, we endow M. K with the L 1 -norm: |x| = Y^k=i \ x k\- Let (x n ), n G N, 
be a sequence in N K such that lim„ \x n \ — oo. For t <E [0, T], we define: 

|a;™| 

To simplify the notation, for f e 1 + , we set W(t) = W(\t]). 

Let B([0,r],R K ) denote the space of cadlag functions from [0, T] to R K endowed 
with the usual Skorokhod topology. This means that the distance between the functions 
/i and /2 is defined by 

<M/i,/2)=inf sup {\g(t)-t\+p(f 1 (g(t)),f 2 (t))}, 
te[o,T] 

where p is the T^-metric in R K and the outer infimum is taken over all monotone 
continuous functions g : [0,T] — > [0, T] such that g(Q) = and g(T) = T. Denote by 
, oo)) the space of _R^-valued cadlag functions on [0, oo) with the metric 

dr(/i,/a) 



T=l 



+ <tr(/i ) / 2 )' 



Note that X n e D([0, T],R K ), for all n. 



Lemma 1. fi) _Fbr any sequence x n such that \x ll \ — > oo, i/ie family {{X n ),n G N} 
almost surely has a compact closure in the Skorokhod topology, and any accumulation 
point z of {X n } is almost surely continuous. 

(ii) Function z is Lipschitz with the constant maxj^^ Xi,K}. 

Proof, (i) One can easily obtain a proof of this statement along the lines of the 
proof of [3], Theorem 4.1 or [T7], Theorem 7.1. Formally, the proofs of the mentioned 
theorems are given for multi-class networks. However, as pointed out in [10] . the 
tightness of such families holds under weaker conditions (see [10], Assumption 2.19). 

(ii) Since Gi are bounded by 1 from above, 

( - l> n |t| ) 

^ max E B k ,K(t-a + l/\x»\)\, 



\x 

1 1 fc^nx^is] 
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where Bk is the total number of arrivals at time k. Sequence {_Bfc}fc g N consists of i.i.d. 
random variables with F,B^ = KA. By the law of large numbers, the result now follows 
if we let n — > oo. 

Definition 1. Any accumulation (in the Skorokhod topology) point z of the sequence 
X n is a fluid limit. The collection of all fluid limits is called the fluid model. 

It follows from the definition of X n and z that |z(0)| = 1 and that Zi{t) > for all 
i = 1, . . . , K and for all i. 

Corollary 1. The sample-path trajectories of fluid limits are self-similar. More pre- 
cisely, for any fluid limit z and for any u > such that P(|z(u)| > 0) > 0, the random 
process {z(t),t > 0} with conditional distribution 



is also a fluid limit on the event \z(u)\ > 0. 

This result may be obtained along the lines of the proof of Stolyar [17], Lemma 6.1. 
However, the remark given in the proof of Lemma [1] (i) also applies here. 

Definition 2. We say that the fluid model is stable if there exists a deterministic time 
to and e > 0, such that \z(t)\ < e a.s. for t > to, for all fluid limits z. Due to the self- 
similarity of fluid limits, this is equivalent to saying that there exists a deterministic 
time t such that \z(t) \ = a.s., for all fluid limits z. 

Definition 2 of fluid stability has become standard and appears in most papers 
dealing with the fluid approximation method. 

2.2. Fluid model criterion for stability 

In this subsection we formulate a stability theorem for fluid limits which will imply 
Theorem Q] (i). 

Lemma 2. // the fluid model is stable then W is positive recurrent. 

Proof. One can again obtain a proof of this assertion by following the lines of the 
proofs of Dai [3] or Stolyar [17] which are given for multi-class networks. 

By Lemma [2] Theorem [T] (i) will follow from the next statement. 




A Random Multiple Access Protocol 



Theorem 2. If X < e 1 /V , then the fluid model is stable. 
Our proof of Theorem is given in Section 2] 

2.3. Fluid limit evolution equation 

In what follows we write 



tpi(t) =&(*(*)) 



Zi{t) 



E *,•(*)' 



(3) 



Theorem 3. Tafce any /Zmd /jmi£ z and fix t > 0. Assume that Ejey z j W > / or 
i. T/iera is differentiable at point t ift>0, and has a right derivative at point 
t ift = 0. Moreover, 



- £ Pi(t) 

z j (t) = A i - ¥ j i (t)e = Ai - Gi(js,(t)), 



(4) 



where Zi(t) is the (right) derivative. 



Under the assumptions of the Theorem, this differential equation admits a unique 
solution. 

Fluid limits with an initial condition on the boundary (Ejev; z j(®) ~ f° r some *) 
are analysed in Section [3] 

Proof of Theorem^ Two cases are possible: either Zi(t) > 0, or zi(t) — and 
J2j£Vi z i(t) > u< We treat these two cases separately. 

(i) Suppose first that Zi(t) > 0. To treat this case, we need the following technical 
result. 

Lemma 3. There exists C > such that \Gi(x) — Gi(x)\ < min(l, C/xi) if Xi > 2. 
Proof of Lemma\^ Using that \e~ Vl — e V2 \ < \yi — y%\ for all J/1,2/2 > 0, we obtain 

1 



\Gi{x) - Gi(x)\ < 



In 1 - 



E 

iev; 



1 



For every j, denote j/j 



From the inequality | ln(l — y) + y\ < 



(5) 



2(1 -2/) 5 
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for y £ (0, 1) we obtain that 

\Gi{x) -Gi(x)\ < yi + 



2(1 -y.f ' S 2(1-^)2 • 



The statement of the lemma now follows from the inequalities 
Vj < 1/aii, ajj-J/j < 1 and yj < 1/2 

for all j €Vi. 

Assume now that t = (the result for an arbitrary t follows from the self-similarity 
of fluid limits). Let x = 2(0). Suppose that s < Xi. Let k < \x n \s, then Wf (k) > 
x™ — k > \x n \(x™ /\x n \ — s). Hence, Wf (k) > 2 for k < \x n \s, for large enough n. 

We need to show that lim s ^ — — ^ = Aj — Gi(z(0)). Consider the decompo- 



sition 



l - ]T (wf(k + i)-wf(k)) 



k=0 
\x"\s]— 1 



fc = 



Wf{k + l)-Wf{k)\W xn {k) 



£ (w? B (A + l)-E[Wf n (* + l)|Wf B (A)] 

|»J-i l>"kJ 
— £ (A,-G i (^"(fc)))+-- ^ ^ (6) 



k=0 



k=l 



where 



Dl = Wf (k) - E [Wf (k) | W x " (k - 1) = A< (k) - Xi + qi (k) - E ( qi (k) \W* n (k - 1) 



with % (fc) = 7(^-1) = l)I\ ]eVMl] I(N J (k-l)=0). Wehave — Ei':i' SJ (^(fc)- 
Xi) — > a.s. as n — > oo. So we can apply Theorem VII. 3 of Feller [8] (with bu = 1/fc) 



to deduce that 



±- £ (ft(*)-E( gi (fc)|W* n (A-l)))-fO a.s. 



(7) 



fc=i 



as n — > oo. 

It remains to find the limit of the first term in the RHS of equation © . Decompose 
this term as follows: 
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jja^lsj— 1 Ma? 71 |s_| — 1 

jk\ g (^-^"(*))) = ]^ E (^-^(^(^Jjj+^n), 
where by Lemma [3] 

1 1 k=0 11 fe=0 1 

as n — y oo uniformly in s < Xi. Further, from the uniform convergence of X n to z and 
the continuity of G we deduce that 

l>"W-i „ / ( k \\ 
*">- "<°>= A, -lin, £ *VW>. 

S n-^oo \x n \S 

Since Y]\j— 1 Gj(z(fen-)) is a Riemann sum of a continuous bounded function, it 
converges to J* Gi(z(u))du, so 

lim «(')-«(°> = Ai - lim /°^(^")) du = At - g 4 (*(0)). (8) 

s^O s s^O S 

(ii) Now consider the second case, Zi(0) — and J2jeVi z i(®) > 0. Notice that 
Gi(zi(0)) = 0. In view of equations © and (JT]) it suffices to show that 

Jx"|«J— l 

s^0+ n.— «x> \x n \s 

By the assumption, there exists j G such that Zj(0) — liru n _voo x™ /\x n \ > a > 0. 
Let e > 0, then there exists no such that for all n > no, Xj/\x n \ > a and < e. 

Fix < s < a and e < a. Then for n large enough, W] 1 (fc) < e\x n \ + Ai(0, k), 
Wf{k) > a\x n \ - k and 

Wf{k) 



lim lim JL £ Gi(W n (k))=Q. (9) 



Gi(W* (fc)) < 



< 



Wf{k)+Wf{k) 
e\x n \ + Ai(Q,k) 



(a + e)\x n \ - k 

By the strong law of large numbers, limt_>+oo Ai(0,t)/t — A; a.s. Let A > Xi. We 
may choose \x\ and n large enough and then fco such that for fco < k < s\x n \, 

e\x n \ + \k 



G t (W x (A;)) < 



(a + e)\x n \ - fc' 
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and 



1 1 fe=0 1 



- £ — S 



|»J- 

E 

fc=0 



a|ic n | — fc 



Direct computations show that 

1 L|;r " |sJ_1 fcA s 
lim : r > — : - = — A(s + aln(l )). 

n^oo \x n \ a\x n \-k y y a" 

Then 

|»J-i 

lim „ su P^[ E ^(W" (*))<- T i_-A( a + a ln(l -1)) a.s. 
Since the last inequality holds for all e > and A > Aj, we have 

|x"|sj-l 



limsup^- V Gi(W^"(fc))<-Ai(a + aln(l--)). 

n |x | k=Q a 

It then follows immediately that 

Ll^isj-i 

lim lim sup 1 — r V Gi(W xn (fe)) = 0. 
s^0+ n \x n \s ^ 

The proof of Theorem 3 is now complete. 

We also need a further result that may be deduced from Theorem [3l Let H — {x G 
R-'*' : Xj > for all i = 1, . . . , if} be the interior of the positive orthant. 

Lemma 4. Assume that z(0) G if, i/iert either 

(i) there exists c such that z(c) — and z(t) G -ff /or aZZ i G (0, c) or 
z(t) remains in H for all t > 0. 

Proof of Lemma^ Restricted on the open set H, the RHS of equation ^ is a 
continuous function. Therefore, the solutions of equation @ are locally uniquely 
defined as long as z(t) remains in H. Now, suppose, to the contrary, that t i— ► z(t) 
leaves H at time c at a point y — limt_> c _ z(t) G dH\{0}. 

Let 0^ = limsup t ^ c - <j>(z{t)), ai G [0,1]. Since y ^ 0, there exist «i and 12 such 
that yi t = and y^ 2 > 0. The connectivity of G implies that there exists k such 
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that yk = and X^gVs, Uk > ® (for that, consider any path from i\ to iz). Hence, 
cik = and lim t _> c - Fk{<p(z(t))) = A > 0, this implies that 1 1— > Xk(i) increases in a left 
neighbourhood of c, which contradicts = lim t ^ c _ z(i) = 0. 

Lemma U] implies that for an initial condition in H the fluid limit z(t) remains in H 
or finally reaches at time c. It also implies that if z(0) = lim„ x n /\x n \ £ H then the 
fluid limit is deterministic. 



3. Properties of the fluid limit on the boundary 

In this section we work with the general case (we no longer make the symmetry 
assumptions). 

Conjecture 1. We conjecture that all coordinates of any fluid limit z have right 
derivatives at point (even if there exists i such that Xj — Zj(0) — for all j £ Vi). 
We also conjecture that the right derivative z (0) of a fluid limit {z(t),t > 0} does not 
depend on the sequence x n and only depends on x = lim„ x n /\x n \. If this is true, then 
all fluid limits are deterministic functions. 

In this Section, we prove a weaker statement which implies that the boundary of 
the positive orthant does not play any role in determining the stability conditions for 
the fluid model. Denote 

Th = mf{f > : \z(t)\ < h}. 

Denote also A* = min{Ai, . . . , A^} > 0. Since |z(0)| = 1, max.;Zi(0) > 1/K. The 
inequality z i (t) > A* — 1 for all i and t also implies that 

^ * W^xry (10) 

Theorem 4. There exist positive constants b and Sq such that, for any e £ (0,£o), 
minj Zi{t) > be for any t £ [ce, Ti_ e ) where c = 1/K(1 — A*). 

From the self-similar property of fluid limits, the following corollary is immediate. 



Corollary 2. For any h > 0, z,(t) > for all < t < Th and all i. 
Lemma [4] and Corollary [2] imply the following. 
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Corollary 3. Assume that \z(0)\ > 0, then either 

(i) there exists c such that z(c) = and z(t) remains in H for all t £ (0, c) or 

(ii) z(t) remains in H for all t > 0. 

We do not present proofs of these statements here as they are rather obvious. 
The rest of this Section is devoted to the proof of Theorem 3J We begin with two 
technical lemmas. 

Lemma 5. There exist positive constants K\ > 1 and Ki such that, for any fluid limit 
z, if Zi(t) > K\Zj{t) for two neighbouring nodes i and j , then Zj(t) > Ki. 

Proof of Lemma [3J The existence of Zj (t) is guaranteed by Theorem^ Since Efcev- z k (*) > 
Zi(t) > 0, we have 



z'( t ) > A* - - Zj{t) , - > A* — - , - >K- 



and we may take K\ = 2/A* — 1 and Ki = A*/2. 

Lemma 6. There exist constants C\ > C2 > such that for any h > one can choose 
hi > such that if \z(0)\ > hi and mirii Zi(0) > Cih then min,z,(t) > Cih\ for all 

t<T h . 

Proof of Lemma Denote by D the maximum graph distance in Q (the diameter 
1 C 

of Q). Put Ci = fttt and Co — — ft^t- We may prove Lemma [5] for h = 1. The 

result for an arbitrary /i follows from the self-similarity of fluid limits. 

It is sufficient to show that for any t < t\ if minjZj(t) > Ci then there exists 
< s < 00 such that 

minzi(t + s) > Ci (11) 

i 

and 

minzi(w) > C2 for all t < u < t + s. (12) 

i 

Indeed, assume that pi p -(|12 p hold and Lemma [H] is not valid. Then there exists t < t\ 
such that minj z i(t) < C2. It then follows from the continuity of fluid limits that there 
is a last moment v < t when miiii Zi(v) > C\. However, it follows from (JTT|)-(|T2J) 
that there exists s > such that mhii Zi(v + s) > C\ and miiii Zi(u) > Ci for all 
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v < u < v + s. Clearly, v + s < t, which contradicts our assumption that v is a last 
moment before t when min^ Zi(v) > G\. 

Now let t be any time such that t < t\ and minj Zi(t) > C\. Note that max^ Zi(t) > 
1/N = C\K® +1 since t < t%. To simplify the notation, assume that Zi(t) — max, z,(t). 
Let T be such that z\{t + u) > C\K^ for all < u < T. Again note that z^u) > A* — 1 
for all i and it. This implies that 

T ^ C^K™ Kf) C 1 Kf{K 1 1) 
1 A* 1 A* 

Let d be the maximum distance in Q from node 1. Clearly, d < D. For j = 1, . . . , d, 

denote by Aj the set of nodes at distance j from node 1 . 

We show that there exists < s < T such that (jXTJ) and (TT2]) hold. First, we show 

that minzi(M) > C2 for all t < u < t + T. Note that Zi(u) > C\ for all i G A\ and 

t < u < t + T . Indeed, assume that there exist i € A\ and t < u < t + T such that 

Zi(u) < C\. Then, by continuity, there is a last moment t < u± < u such that Zi(u\) > 

C\. Lemma [5] implies that z i (u\) > K2 > and hence, there exists time 1*2 > U\ such 

that Zi(u-i) > C\, but this contradicts our assumption on u\. Using induction and 

following the same arguments, we can show that Zi{u) > C\jK\ 1 for all i G Aj and 

t < u < t + T for any j = 1, . . . , d. Hence, mi ni Zl {u) > CxjK^ 1 > d/K^ 1 = C 2 

for alU < u < t + T. 

Let us now show that there exists < s < T such that pip holds. For every 

j = 1, . . . , d, denote by tj the time needed to achieve the level CiK-^ 1 starting from 

C (K^~^ — 1) 

the level C\j ' K[~ l and moving with speed K2. Clearly, tj = 1 — . Note that 

d d 

PT|) and (JT3J) hold with s — ^ tj if T > tj- Indeed, minjg^ Zj will achieve the 

j=l 3=1 

level Ciiff -1 not later than at time t + ti and will not become smaller than this level 
before time t + T, since all nodes in Ai are neighbours of node 1 and Zi(u) > for 
all t < u < t + T. Note also that mhijg^ Zj will become greater than C\Kf~ 2 not 
later than at time t + 1\ + ti since it cannot become smaller than C\ j K\ before time 
t + ti, and after this time it is either greater than CiK^~ 2 or grows with a speed of at 
least K 2 (this follows from Lemma [5] and the fact that any node in A 2 has a neighbour 
in A\). We can continue these arguments to prove that rrrhx/g^ Zj will become greater 

d d 

than C\ not later than at time t + U if T > ^ U. 

i=l i=l 
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Note that 

^ . C^Kf- 1 - 1)(1 + K x + . . . + Kf- 1 ) _ C^K*- 1 - \){Kf - 1) 



that D is the maximum graph distance of Q. Take e such that — 77775 — — < 



ti 1 KiKf- 1 K 2 Kt X {K x - 1) 

< CxjKf-1) < -1) 

" - 1) " # 2 («i - 1) • 1 ' 

If we take X 2 = A*/2 and Xi = 2/A* - 1 then (1 - K)/K 2 = K\ - 1. Note also 

d 

that in this case K x > 2. It now follows from (H]) and HU) that T > £ U- 

i=l 

One can see from the proof of Lemma [5] that the following (stronger) result holds. 

Lemma 7. For any hi > there exists h 2 > suc/i i/iai /or any h 2 < ft 2 i/iere exists 
< /13 < «2 smc/i £/ia£ i/ z(0) > ft-i and miniZi(O) > /i 2 then minj^j(t) > /13 /or all 

t<T hl . 

Remark 1. Lemma [7] is valid with h 2 = — —fttt- 

Proof of Theorem^ The proof of Theorem H] is similar to that of Lemma Recall 

K 2 {K X - l)e 

(KP - 1) 

fttt for all e < en and let a = * ^— . Then ae < ftt-t, and, in view of 

Lemma [7] and Remark [TJ it is enough to prove that min^ Zi(ce) > ae. 

Note that maxiZ^O) > 1/K. Assume that zi(0) = maxiZi(O). Let T be such that 

zi(u) > asKP for all < u < T. Then z'^t) > A* — 1 implies that 

T . VK- aeKP _ 1 - KaeKP 

~ 1-A. " K{l-K) ■ (15) 

Again let d be the maximum distance in Q from node 1. For j — 1, . . . , d, denote 

by Aj the set of nodes at distance j from node 1. For every j = 1, . . . , d, denote by 

tj the time needed to achieve the level aeK^~^ starting from the level and moving 

aeK d ~ } D 

with speed A 2 . Clearly, tj — — . Denote Xi = ^ tj- Note that 

K 2 J= i 

ae(KP - 1) e 
Tl= A 2 (A 1 -1) = A(T-A0= C£ - (16) 

Following the same arguments as in the proof of Lemma [6j we can show that 
mini Zi(cs) = min* Zi{T\) > ae if T\ < T . 



It remains to prove that T\ < T. This follows from (fl5|) . (fl6|) and the inequality 
1 -e 

ae < 



KK{ 



D+i ■ 
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Remark 2. Denote by v(s,h,b) = ini{t > s : \z(t)\ < h or mirij 2j (t) < b} the 
time (after moment s) of the first exit from the set {\z\ > h} fl {mini #j > 6} after time 
s. Theorem [4] implies that there exist b > and s > such that ri_ e = v(se, 1 — e, te) 
for any initial condition z(0). 

4. Proof of Theorem [2] 

In this Section we first present a proof of our main stability result and then formulate 
its generalisation. Recall that Theorem 1 follows from Theorem [2] and Lemma [2] 

Recall also that here we deal with the symmetric case. We need to prove that there 
exists a deterministic time to such that for all fluid limits, z{t) = for t > t$ a.s. 

Lemma 8. // Zi(t) > for all i = 1, . . . , K , then 



and hence, if A < 



v ' 



for some e > 0. 

Proof of Lemma . Clearly, it is sufficient to prove the inequality 

z l (p l exp I - £ ^ I 
» ieVi _ e" 1 



> -FT ( 17 ) 



k 

where we slightly abuse the notation from ^ by writing instead of zi(t) and ipi 
instead of ifi(t). We can write the LHS of the previous inequality in the form 

i 

Zi 1 

where pi = — — , Hi = — X ^Pj ~ m — an d f( z ) — e * ■ Function / is convex 

X z k jeVi fi 

k 

and J^Pi = !. hence Y,Pif(Vi) > f(52PiVi) and 



V A; fc 
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Now consider ^2 — X fj an d S v^" - ^ n — separately: 

i E z k jeVi i Z_^ z k fi 

k k 

' ?^S W = ^ _ = ^ _ = 5^ = 

(we use the identity ^ = Zj and the symmetry of the neighbourhood relation: 

j e Vi iff z 6 V$.) 

Note that the logarithmic function is concave. Hence 

£ * in 1 < ta *lL m f = m \^-\ = mv. (20) 

Inequality JTTJ) follows now from (JTSJ), (JTSJ) and (|20]). 

Proof of Theorem^. Corollary |2] implies that Zj(t) > for all z = 1, . . . , K and all 
< t < inf {u : z(u) = 0}. Then we can use Lemma [SJ Note also that ^ x% > 4 : 

for any positive values of {xi}. Hence, Lemma [5] implies that 



Then 



and the result follows. 



/£><(*) <-e/2, 



Remark 3. By applying methods of [TH], we can get a similar (but less explicit) 
stability result in a more general situation. Assume now that the system may be 
asymmetric, i.e. that values of A.; may differ for different i and the graph Q may be 
irregular. 
Let 

M = {[i : ^ =p;e~^ ev » P3 ,i = 1, . . . ,K, for some p = (p x , . . . ,p K ) with p t >0}. 

One can show that the vector (tp±, . , . , (px) with ipi = — — maximises the 

EjeVj z j 

K 

function Zi In /ij over all vectors /j, S M . Based on that, one can obtain the following. 

i=l 
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Theorem 5. Assume that there exists a vector /i 6 M such that A < fi component- 
wise. Then the Markov chain W n is positive recurrent. 

A proof of Theorem [5] follows the lines of the proof of Theorem 4 in [18] . 

5. Rate of convergence 

In this section, we again consider the symmetric case. We will obtain power rates 
of convergence of W n to its stationary distribution in the total variation norm. We 
expect that one can similarly prove the geometric ergodicity of the underlying Markov 
chain given the light-tailedness of the distribution of the increments {A(n)}. 

Define the total variation distance between distributions ~k\ and 7T2 by 

IM-) -7r 2 (-)|| = sup 
lffl<i 

Theorem 6. Assume that A < e^ 1 /V and ~EiAi{n) p+1 < oo for some p > 1 and for 
all i = 1, . . . , K and n. Assume also that P(Al(1) = 0, ^(1) = 0, Ak(1) = 0) > 0. 
Then 

lim n p ||P™(izv) -tt(-)|| = 0, xeN K , 

n — >oo 

where P"(x, •) is the distribution ofW x (n) and ir(-) is the stationary distribution of 
W. 

Proof of Theorem The proof of Theorem[S]is based on the following lemma which 
is an analogue of Proposition 5.3 of Dai and Meyn 0]. 

Lemma 9. Assume that the conditions of Theorem [6] are satisfied. Then, for some 
constants c < oo, 5 > and a finite set C, 

(ra{S) \ 

for any x £ N K , where tc(S) — mm(n > S : W(n) £ C). 

Proof of Lemma\^ The proof of Lemma [9] follows the lines of the proof of Proposi- 
tion 5.3 of @]. 

It follows from Theorem [2] that there exists to such that 

lim W 'flf°> = a.s. 

|x|->oo |x| 



g(y)ni(dy) - / g(y)Tr 2 (dy) 
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f |W x (b|£o)r° +1 ] 

Note also that the family of random variables ^ ■ — j^-j )■ is uniformly inte- 



grable, since 

\w*(\ x \t )f +1 < (E^° EL^m) p+1 < tp+1 Elt° (EL^M 

\x\t 



• |P+i 



and the family 



p+i 



Em=°0 (E^lAjj 

|a;|to 



p+i 



is uniformly integrable. This is guaran- 



teed by the existence of ~EAi(m) p+1 for alH = 1, . . . , K and for all m. Hence, 

lim V[\W*(\x\to)\r^) q 



I p+i 



Choose L such that 



Efl^dxlto)!^ 1 ] <\\x\ p+1 



(21) 



for > L. Define, as in the proof of Proposition 5.3 of [4], the sequence of stopping 
times do = 0,<7i = t(x), and dk+i = &k + 8 ak ai,k > 1, where t(x) = tomax(L, \x\) 
and 8 is shift operator on the sample space. We assume that to is an integer. The 
stochastic process Wk = W{ak) is a Markov chain with the transition kernel 

P(x,A) = P(W x (t(x)) e A). 

Now (|2"Tj) implies that 

E {\W^ +1 - \W \p +1 \W = x} < -^\x\p +1 +bl c (x), 

for the set C = {x : \x\ < L} and for some constant b. The Comparison Theorem 
(Meyn and Tweedie [13] . p. 337) yields that 



E 



fc„-i 

E i^Kor 1 

n=0 



E 



E \wm p+1 

,n=0 



<2{\x\ p+1 +bl c (x)} (22) 



where = min(fc > 1 : W{k) G C}. To prove Lemma |9l we first show that for some 
constant cq 



E 



E \W x (n)\*\T ak 



< c W x (a k ) 



x (rr, \P+ l 



(23) 



which by the strong Markov property amounts to 

t(x) 

Ej2\W x (n)\ p < c \x\ p+1 . 

n=0 
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This follows from the fact that 



t(x) t(x) I n K \P t{x) I t(x) K 



n=0 



n—0 \m— i—1 



n—0 \ m—0 i—1 



a.s., and from our assumption that EA^to) < oo for all i = 1, . . . , K and for all m. 
Substituting (|23[) into (j2"2"|) , we have 

~<?k + l 



E 



,fc=0 



E i r wr^ 



<c|.t| p+1 . 



By Fubini theorem and the smoothing property of the conditional expectation, the 
LHS is precisely E [X)n=o(l + l^ a: ( n )l p )] ■ The proposition now follows from the fact 
that a kt > T C {t L). 

We now apply Proposition 5.4 of ;4 with t = 1. In our case, it gives the following 
bound: 

E {V(W(l)) - V(W(0))\W(0) = x} < -f(x) + k (24) 



with V(x) = E (j2l=o ] \W x (n)\P^j and f(x) 



Further, Lemma M implies that V(x) < c\x\ p+1 . Now ([M]) yields that 

E{V(W(1)) - V(W(0))\W(0) = x}< V(x)^ + bl c 

for the set C = {x : \x\ < L} and for some constant b. The result now follows from 
Theorem 2.5 of Douc et al. [5]. 



6. Local stability of fluid limits on the positive orthant 

In this Section we investigate the behaviour of the solution to the differential equa- 
tion (j4|). In particular, we show that if the input rate A is sufficiently small, then the 
diagonal is locally unstable. 

6.1. Orbits of the fluid limits 

Recall that H = {x e R K : Xi > for all i = 1, . . . , K} is the interior of the 
positive orthant and 1 = (1, 1). For z(t) in H, the differential equation ((4]) may be 
restated in a closed form as 

z'(t)=F(cb(z(t))), (25) 
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with Fi(x) = X — x i e^^-' jev i Xl . Let A = {x G H : x\ — x<i... — %} be the main 
diagonal and C u = {x G H : \x/\x\ — ^/K\ < u}, u > 0, C u a cone with direction A. 
We note that the main diagonal is an orbit of the equation F((f)(d)) = (A — e _1 /V)l. 

Let A be the adjacency matrix of £/ and {y\, its eigenvalues with Vi < Vi+i. 
The spectral gap 7 is defined by: 

7 = min(^ - Vi) = v K - v K ~i- 

Note that since Q is (V — l)-regular, vk = V. The main result of this section is the 
following. 

- 1 2 

Theorem 7. If X > j tt^(1 — y^)i then there exists u > swc/i £/ia£, for all solutions 
t — > 21(f) 0/ equation 125)) with the initial condition in C u , 

lim 0(z(t)) = 
£ — >+oo 

—1 2 

If X < ^y^(l — ys); fien f/ie diagonal is locally unstable. 

Theorem [7] will be proved in the next Subsection. 

_i 2 

Corollary 4. Assume that X > ^rr-(l — 7^) and that z(t) is a solution of 125]) . There 
exists u > such that if z(0) G C U; then 

(i) if X < e~ l /V , then z(c) — for some c > 0. 
(ii) ifX> e-yV, fften z(f) - (A - e" 1 /^*. 

Proof of Corollary^ Let z(t) be the maximal solution with given initial value z(0) G 
H. From ThcoremO lim0(z(t)) = 1/V. Since F is continuous in a neighbourhood of 
l/V, lim t ^ +00 z'(f) = (A - e _1 /y)l. If A ^ e^/V, then the latter implies that, as f 
tends to infinity, 

z(t) ~ (A - e _1 /V)tl. (26) 

Then the second statement of the corollary follows. Suppose now that A— e~ x jV < 0. 
Then, from (|26[) . 2(f) leaves -ff in finite time. Lemma [4] implies in turn that there exists 
c > such that z(c) = 0. So, the first assertion of Corollary [5] is also proved. 

6.2. Proof of Theorem [7] 

The proof of Theorem [7] is an application of the stability theory of differential 
equations. It will be given in the series of technical lemmas. 
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6.2.1. Spectral analysis. We need to consider the eigenvalues of D(F o <fi)(x) for x G A, 
where Df(x) is the differential of / at x. Here F o <$> is homogeneous of order 0: for all 
c> 0, F((f>(cx)) = F((j){x)). Hence, 

D(F o 0)(cl) = c^DiF o </>)(!). 

Since A is an orbit of equation ([23)) . 1 is an eigenvector of _D(i 7 'o^))(l) associated with 
the eigenvalue 0. 

Lemma 10. The eigenvalues of D(Fo(f>)(\) are (0, 771, •• ■ , 77.R--1) withrji — ~^r(V — 
VK-i) 2 ■ In particular, rji < /or alii > 1. 

Proof of Lemma \TUl A direct computation leads to: 

y2 II J — I, 

llMFo. >){!)),; = \ ^.\ViUVj\ if jeVi\{i}, 

-&\VinVj\ if jiVt. 

Then £>(F o </>)(!). 1 = 0. To show this, let M = -eV 3 D(F o 0)(1). Using the 
equality \Vi LI Vj \ = 2V - \Vi n Vj\, we deduce that: 

(mi), = v(v - 1) - 2F(y - 1) + 53 \Vi n I = ^ 1*5 n ^1 - ^ 2 = °- 

Let -E denote the identity matrix and j4 the adjacency matrix of Q. Since (A )ij — 
\Vi PI Vj\, we have the following decomposition: 

M = V 2 £ - 2VA + A 2 = (A - VE) 2 . 

The matrix A is irreducible since Q is connected. Thus (A — VE) is an ML-Matrix 
(refer to Seneta |16j). In the graph theory, this matrix is called the Laplacian matrix 
of Q. From Corollary 1 of Theorem 1 in Seneta [16], the spectral radius of A is V. 
Theorem 2.6 (d) of [H] implies that dim K er {A — VE) — 1 and that all non-zero 
eigenvalues of (A — E) 2 are positive reals (recall that the spectrum of A is real). 

6.2.2. Orbit of ip ° z. We define: 

K 

V = {xeH : ^a; i = l} = irn(l,-)- 1 ({l})=V'(^), 

i=l 
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where ^p(x) — x/\x\. S is clearly a C°°-convex manifold of codimension 1. Introduce 
the following differential equation on E: 

y' = D^(y)F(4>(y)) - a(y) (27) 

with an initial condition y(0) in S. Here a is a C°°(S) function and a(y) G T y (Yj) the 
tangent space of £ at y. The next step is to compare the orbits of equations (|27p and 
([23)1 . The next lemma asserts that the orbits of the solution of equation y' = a(y) and 
f/)oi are indeed equal (here 1 1— > z(t) is a solution of equation l|25[)). 

Lemma 11. Let z(0) belony to H and let z(t) be the maximal solution of equation 125\) . 
Let y(t) be the maximal solution of y' = G{y), with the initial condition y(0) = ip(z(0)). 
Then y{t) is defined on R+, and there exists an increasiny continuous bijective function 
fi : K_|_ — > M + such that 

y o fx = ip o z 



Proof of Lemma[TJ[ For any initial condition in H, we have F(tfi(z(t))) < Al. This 
is clearly true if z(t) G H. If z(t) £ H, then z(t) G A n H c , by Lemma |H Thus 
F(z(t)) = (X-l/Ve- 1 )! < Al. It follows that \z{t)\ = < KM+Th-i z j( )- 

Suppose now that z(t) G H for all t. Then the integral (j2f=i Zj{s)j ds 
diverges. By the intermediate value theorem, we deduce that there exists an increasing 
continuous function v such that 

'"(«) ds 

= t for all t > 0. (28) 

In particular, 

K 

i/(t) =$>("(*))• 

3=1 

Let w = %jj o z o with iu(0) = ib(w(0)) = 2/(0). We have 



ds 



s=v{t) 



K 



Y,zMt))\mz(Ht)))F(w(t)). 
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The function ip is homogeneous of order and thus Dip(cz) = c 1 Dip(z), for all 
c> 0. Then 

w\t) = DtP ( / Kt)) ) F(w(t)) 
= G(w(t)). 

The solution of the differential equation is unique, therefore w(t) = y{t). The lemma 
is proved with /j = v . 

6.2.3. Local stability of ip o z. Clearly, yo = 1/K is an equilibrium point of equation 
(|27|) . In the next lemma we prove that this equilibrium is locally stable. 

— 1 2 

Lemma 12. If X > S— (1 — y?), there exists u > such that, for all solutions 1 1— > y(t) 
of equation with \y(0) — yo\ < u, 

lim sup \y(t) - 1/V\ = 0. 

t ^ + °°y(o)es:|y(o)-ao|<u 

Proof of Lemma ll^i We denote by Dck(jj)\t (s) the differential of a at y restricted 
to the K — 1 dimensional subspace Tj,(E). It is known that if all the eigenvalues of 
Da(yo)\T y C£) nave a negative real part, the local stability follows (see, e.g., Coddington 
and Levinson [2]). Let D 2 ip(y)(- 7 •) denote the second differential of ip at y, seen as a 
bilinear mapping. We have 

Da(y) = D 2 ^(y)(F(cf>(y)), •) + D^(y)D(F o 0)(y). (29) 

The first term in (I29[) is a matrix and its entry (i,j) is equal to 

Clearly, F(<f)(yo)) = (A — e _1 /V)l. Then a straightforward computation gives 

D 2 i/j(y )(F^(y )), •) = (A - e- x /y)(J - 

where i5 is the identity matrix and J is the matrix with all its entries equal to 1. We 
also have Dtp(yo) = E — J/K. Finally, equation ([29)1 can be rewritten as 

Da(y ) = (E - J/K) ( D(F o 0)(y o ) - (A - e' 1 /V)E ) . 
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The matrix (E—J/K) commutes with all symmetric matrices and has two eigenvalues: 
eigenvalue 1 (with multiplicity K — 1) and eigenvalue (with multiplicity 1, associated 
to the eigenvector 1). By Lemma [TOl the eigenvalues of D(F o 4>)(yo) — (A — e~ x /V)E 
are 

(M = -e -1 (V - v K -i) 2 /V 3 - X + e^/V for 0<i<K-l. 

The eigenvector associated to /io = A— e^ 1 /V is 1. Thus we have proved that A — e _1 /V 
is an eigenvalue of multiplicity 1 for Da(yo) and that the other eigenvalues are (ni)i>i- 
These eigenvalues have negative real parts if and only if \l\ = — e^ 1 ^ 2 /V 3 — A + 
e _1 /V < 0, which is equivalent to A > e _1 (l — r y 2 /V 2 )/V. The vector space generated 
by the associated eigenvectors is precisely the tangent hyperplane T 2/0 (S) = I' 1 , the 
hyperplane orthogonal to 1. 

Now we can prove Theorem [7] Let \z(0)\ e C u and y(0) = z(0)/|z(0)|. Then, by 
Lemmas Qj] and [T^l 

lim ^(z(t))= lim y(»(t))=l/K. 

t — »+oo t — >+oo 

In particular, <fi(z(t)) tends towards \/V as t tends to infinity. 



7. Absence of attraction to the diagonal in one particular case 

It has already been pointed out in Section 1.3 and in the previous Section that if A 
is too small, the main diagonal may not be locally stable. In this Section, we present 
an example of a graph with locally stable sets of parameters which do not belong to 
the main diagonal if A is sufficiently small. 

Consider a graph Q with 4 vertices placed on a circle. Number the vertices 1, 2, 3, 4 
clockwise and assume that each vertex is linked with its 2 neighbours (so that, for 
example, vertex 1 has links with 2 and 4). In this case, K = 4 and V = 3. 

For the fluid limits associated with this graph, consider equation (|27|) . It can be 
rewritten in the form 

K 

y'i(t) = (A - <pi{t)e- w<*)) _ y. J2 (A - <p k (t)e- E ^ fc wW) , i = l,...,K. 

fe=i 

We are interested in the so-called stable points of the latter system of differential 
equations, i.e. points for which the RHS's of all the equations above are identically 0. 
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So if y(0) is such a point, y(t) stays at this point for all t > 0. Clearly, one stable point 

e- 1 ( 7 2 \ 

is (l/K, . . . , l/K), which corresponds to the diagonal. However, if A < 1 ^ 

V \ V 1 ) 

(= 5/27e~ 1 in our case), then there exist other stable points. 

Take yi (0) = 7/2 (0) and y 3 (0) = j, 4 (0). Since ^(0) + j/ 2 (0) + y 3 (0) + ^(0) = 1, the 
equality y 3 (0) = (1 — 2y 1 (0))/2 holds, and the system of differential equations at time 
t = reduces to just one (i.e. any) of them. One can show that, for any A < 5/27e _1 , 
the RHS of this equation equals at three different points: at j/ 1 ) = (1/4, 1/4, 1/4, 1/4) 
and at two others, say j/ 2 ) and . One can find approximate values of these points 
numerically. For instance, if A = 0.001, then w (0.01,0.01,0.49,0.49) and y@) w 
(0.49, 0.49, 0.01, 0.01). Numerical results also show that these points are locally stable. 



8. Extensions of the model 

8.1. Random neighbourhood 

In this Subsection we consider a possible extension of our model. Assume there 
is a fixed number of points 1,...,K and a set of undirected graphs {^• 7 } _ 1 each 
having points 1, . . . , K as its vertices. Assume that at each time n the neighbourhood 
relations are given by the graph where r\ n are independent identically distributed 
random variables taking the value j with probability pj . The need to consider such a 
variability of ncigbourhood relations may be justified by, for instance, the fact that a 
change of environment conditions may lead to a change of the radius and/or direction 
of interactions. 

Denote by Vf the neighbourhood of the point i in the graph Qi and by V/ its 
cardinality. We assume again that the system is regular in the sense that EVj n = V 
for all i. 

Following the proof of Theorem [3l one can show that the fluid limits of the model 
described above satisfy the following differential equation 

L - E 

^) = A-]>> fc ^(i)e jev * fc 
fc=i 

where ^{t) are defined in a natural way. Using the same methods as those used in 

the proof of Theorem [2J it can be shown that the system with random neighbourhood 
e" 1 

is stable if A < — — . 
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8.2. Non-regular graphs with space-inhomogeneous input 

Although Remark [3] provides sufficient conditions for stability in this case, the 
conditions are not easy to verify. Here we give some other conditions that are also 
sufficient for the stability of the system. Assume now that E V/' 1 = Vi , and Vi are not 
necessarily equal. Put V — maxVj. Assume also that E£™ = A,, so that the input 

i 

is "space-inhomogeneous". Put A = maxA^. Clearly, all the results concerning fluid 

i 

limits also hold in this case, and it is easy to see that one can prove the following result. 

e- 1 

Theorem 8. The system described above is stable provided A < 
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